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ABSTRACT 

. We discuss the feasibility of the detection of the 21 cm forest in the diffuse IGM with 

1 the radio telescope LOFAR. The optical depth to the 21 cm line has been derived using 

simulations of reionization which include detailed radiative transfer of ionizing pho- 
tons. We find that the spectra from reionization models with similar total comoving 
hydrogen ionizing emissivity but different frequency distribution look remarkably sim- 
ilar. Thus, unless the reionization histories are very different from each other (e.g. a 
predominance of UV vs. x-ray heating) we do not expect to distinguish them by means 
of observations of the 21 cm forest. Because the presence of a strong x-ray background 
k> . would make the detection of 21 cm line absorption impossible, the lack of absorption 

; i " could be used as a probe of the presence/intensity of the x-ray background and the 

C$ \ thermal history of the universe. Along a random line of sight LOFAR could detect a 

global suppression of the spectrum from z <L 12, when the IGM is still mostly neutral 
and cold, in contrast with the more well-defined, albeit broad, absorption features 
visible at lower redshift. Sharp, strong absorption features associated with rare, high 
density pockets of gas could be detected also at z ~ 7 along preferential lines of sight. 
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1 INTRODUCTION 

The hyperfine transition of the ground state of the neutral 
hydrogen atom, 2 Si/2, is the most direct probe of the content 
and distribution of neutral hydrogen in the universe. Recent 
advances in the construction and planning of powerful radio 



interferometers (e.g. LQFAlQ MW/Q, Sk43) are promising the 
opportunity to make direct observations of the high redshift 
intergalactic medium (IGM), at z ^ 6.5. 

The attention dedicated to 21 cm studies in cos- 
mology so far has been mainly focused on 21 cm to- 
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mography, which would ideally provid e a 3D mapping of 



the evolution of n e utral hydrogen (e . g. Madau et al. 19971; 



2006; 



Dattaetal 



Back et al 



20091: 



Tozzi et al l l2000l; Iclardi fe Madaul 120031; iMellema et all 



20081; iLidz et al.ll2008l ; |Santos et al.ll2008l; 
Geil fc Wvithei 120091; Morales fc Wvithel 



20101 ; I Santos et al.ll2010f ). Although this is an extremely ex- 
citing prospect, it suffers from several severe difficulties (e.g. 
Shaver et ahll 19991; iDiMatteo et af1l2002l. |200 j; Ideser et all 



20081; iGeil et al 1120081; I Jelic et al.ll200S); iBernardi et al.ll2009l : 



Bowman et al]|2009l ; iLiu et alj 120091; IBernardi etal]|2010h . 

the most relevant being foreground and ionospheric contami- 
nation, terrestrial interference, low spatial resolution and the 
requirement that the IGM spin temperature, Ts, is decou- 
pled from the temperature of the cosmic microwave back- 
ground (CMB), Tcmb. 

A valid alternative to tomography is looking at the 
21 cm lines generated in absorption against high-z ra- 
dio loud sources by the neutral IGM and intervening 
colla psed structures, i.e. to search for the 21 cm forest 
(e.g.jCarilli et al.ll2002l; IFurlanetto fc Loebll2002l iFurlanettc 



2006 



Carilli et al 



2007 



Xu et al l 120091 ; iMack fc Wvithe 



20 111 ; lMeiksinll201ll ; IXu et al.ll201ll ). Analogous to the ex- 



tensively studied case of the Lya forest (for a review see 
lMeiksinll2009l ), the 21 cm forest signal can be detected in 
the spectra of high- z radio sources and results from the ab- 
sorption produced by the intervening neutral hydrogen along 
the line of sight (LOS), which removes from the spectrum 
the continuous radiation locally redshifted to the resonance 
frequency. 

Despite the great challenge posed by the requirement of 
suffici ently bright, and h e nce extremely rare, target sources 
(e.g ICarilli et all |2002| ; IXu et all 120091 ; IMack fc Wvithel 
l201ll l the 21 cm forest is particularly appealing as it nat- 
urally bypasses the main limitations expected for 21 cm 
tomographic measurements. Using 21 cm forest observa- 
tions provides several benefits: (i) the IGM is visible even if 
T s = Tcmb, because the CMB does not act as background 
source; (it) the IGM can be spatially resolved on small scales, 
of the order of a few tens of kpc or even lower; (Hi) spec- 
troscopy of bright point sources is technically easier than 
2D image processing, and all proble ms associated with fore- 
grounds and interference vanish (see lFurlanetto| [2006. for a 
more detailed discussion). A further advantage is the close 
similarity with the Lya forest, from which one can adopt the 
well established and robust analysis techniques described in 
the literature. 

In this paper we discuss the possibility of detecting the 
21 cm signal with the radio telescope L0FAR, b y making use 
of the simulations of reionization presented in lCiardi et all 
|2012l . C2012). Differently from previous studies mentioned 
above, the simulations used here are run using a combina- 
tion of hydrodynamic simulations and the detailed radiative 
transfer of ionizing photons through a gas composed of hy- 
drogen and helium. The latter component is particularly im- 
portant for a correct determination of the gas temperature. 
In addition, the expected observed spectrum is calculated 
based on the actual LOFAR configuration. In Section 2 we 
briefly review the basic physics relevant for the 21 cm for- 
est signal. In Section 3 we discuss the adopted simulations, 
while in Section 4 the method for producing mock spectra 
of high-z radio loud sources, as well as the dependence on a 
number of quantities, is presented. In Section 5 we discuss 



the observational feasibility and in Section 6 we give our 
conclusions. 



2 THE 21CM FOREST 

As the physics behind the emission/absorption of the 21 cm 
line has been extensively discusse d in the literature by sev- 
eral a uthors (for a review see e.g. IFurlanetto. Oh. fc B riggs 
2006), here we just write the relevant equations without 
derivation. The 21 cm line corresponds to the hyperfine, 
spin-flip transition of the ground state of the neutral hy- 
drogen atom, 2 Si/2, and to a rest frame frequency V2icm 
= 1.42 GHz. A photon traveling through a patch of neu- 
tral hydrogen at redshift z will see the o ptical depth (e.g. 
iMadau et allll997l ; IFurlanetto et al.ll2006t ): 
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where hh is the H number density, 5 is the gas density con- 
trast, xm is the mean neutral hydrogen fraction, Tcmb is the 
CMB temperature, T s is the gas spin temperature (which 
quantifies the relative population of the two levels of the 
2 S 1 /2 transition), A2i cm = 2.85 x 10~ 15 s _1 is the Einstein 
coefficient of the transition, H is the Hubble parameter and 
dv\\/dr\\ is the gradient of the proper velocity along the LOS 
(in km s _1 ), which takes into account also the contribution 
of the gas peculiar velocity. The other symbols appearing 
in the equation above have the standard meaning adopted 
in the literature. It should be noted that here we do not 
perform a convolution with a Doppler line profile. Although 
we do not expect this approximation to affect our estimates 
of the detectability of absorption lines with LOFAR, we are 
planning to investigate this issue in more detail in the future. 

The photons emitted by a radio loud source at redshift 
z s with frequencies v > z^icm, will be removed from the 
source spectrum with a probability (1 — e _T21cm ), absorbed 
by the neutral hydrogen present along the LOS at redshift 
z = (t y 2icm/i y )(l + z s ) — 1. Analogously to the case of the 
Lya forest, this could result in an average suppression of 
the source flux (produced by diffuse neutral hydrogen), as 
well as in a series of isolated absorption lines (produced by 
overdense clumps of neutral hydrogen). 

Due to the several orders of magnitude difference in the 
corresponding Einstein coefficients, the physical regimes in 
which the 21 cm forest and the Lya forest are relevant, are 
significantly different. As a consequence of the extremely 
low value of A2i cm , the 21 cm absorption is not saturated at 
high redshifts, as is the case for the Lya forest. This makes 
the 21 cm forest signal suitable for studying the IGM in the 
epoch prior to reionization, complementarily to the Lya for- 
est which instead is more appropriate for post-reionization 
studies. Nevertheless, the 21 cm forest detection is in prac- 
tice extremely hard, mainly due to the tremendously weak 
21 cm magnetic hyperfine transition. As can be easily in- 
ferred from eq. [1] the associated optical depth is above 
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0.1, and could therefore produce an observable signal, only 
under severe physical conditions (i.e. highly neutral, dense 
and cold gas). Considering furthermore the paucity of radio 
loud sources at the high redshifts of interest, and their anti- 
correlation with the optical depth due to the fact that these 
kinds of objects are thought to be associated with large ion- 
izing fluxes, it is immediately clear that the detection of a 
21 cm forest signal is an extremely challenging experiment. 



3 SIMULATIONS OF REIONIZATION 

The modelling of mock spectra has been done in two main 
steps. As can be seen from eq. [T] the optical depth and 
thus the absorption along the line of sight to a high red- 
shift source, depends significantly on the gas spin temper- 
ature, as well as on the neutral hydrogen fraction at each 
point along the LOS. The evolution of both these quan- 
tities and their spatial fluctuations are determined by the 
complex process of reionization, which needs to be properly 
modelled. The first step in our spectra modelling pipeline, 
which is described below, is then to simulate realistic reion- 
ization scenarios. The second step is to use the outputs of 
such simulations to extract random LOS and build mock 
absorption spectra of sources at different redshifts. More de- 
tails are given below in Section 0] Unless stated otherwise, 
we assume that the spin temperature is coupled to the gas 
temperature, T. 

The simulations we use are described in C2012. Here we 
just give some basic information and we refer the reader to 
the original paper for more details. 

The simulations are based on a combination of hydro- 
dynamic simulations performed using the parallel Tree-SPH 
code GADGET-3 (which is an updat ed version of t he pub- 
licly available code GADGET-2; see ISpringej|2005h and of 
3D radiative transfer of both H and He fo l lowed by the 
Monte Carlo code CRASH (ICiardi et all I2OO1I; iMaselli et all 
120031 . 120091 ; IPierleoni et~alll2009l ; IPartl et al.H201ll ; Graziani 
et al. in prep.). The hydrodynamic simulations were run in 
boxes of comoving size 35.12/i _1 Mpc, with cosmological pa- 
rameters Q,k = 0.74, Q m = 0.26, Cl b = 0.024ft" 2 , h = 0.72, 
n s = 0.95 and ug = 0.85, where the symbols have the usual 
meanings. A total of 2 x 512 3 dark matter and gas parti- 
cles were followed in the simulation, yielding a mass per gas 
particle of 4.15 x 10 6 /i _1 Mq. After the simulation outputs 
were obtained, the gas densities, peculiar velocities, temper- 
atures and halo masses were interpolated on to a uniform 



128 grid for the radiative transfer calculations. The emis- 
sivity and spectra of the sources were chosen based on a 
semi-analytic model which satisfies a number of observa- 
tional constraints, among them the Thomson scattering op- 
tical depth and the H 1 photoionization rate measured from 
the Lya forest. The spectra adopted are always power-laws, 
with index a=l, 1.8 or 3. Our reference run (£1.2-a3 in the 
original paper) has a = 3. In Figure [1] maps of the tem- 
perature (left panels) and neutral fraction (right panels) are 
shown at redshifts z = 14, 9 and 7, which correspond to 
volume-averaged Hn fractions xhii (gas temperature T) of 
0.029, 0.481 and 0.852 (488, 6020 and 10643 K), respectively. 
This simulation mimics a late reionization model, while the 
other three simulations have been run with a higher comov- 
ing hydrogen ionizing emissivity (see eq. 3 of C2012) and 
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Figure 1. Left panel: maps of temperature at redshift z = 14 (up- 
per panel), 9 (central panel) and 7 (lower panel) for the simulation 
£1.2-a3. Each map represents the central slice of the simulation 
box. Right panel: as the left panel but for the HI fraction. 



reach full reionization at an earlier time. These all have a 
similar emissivity assuming that 30% (70%) of the sources 
have a = 1 (3) (model £ 1.2-al-3), or that all sources have 
a = 3 (model £1.6-a3) or a = 1.8 (model £1.2-al.8). The 
characteristics of the models are summarized in Table 1 of 
C2012. As mentioned, with the exception of £1.2-a3 all the 
simulations have similar comoving hydrogen ionizing emis- 
sivities, so any differences in the reionization histories are 
only due to the different spectral energy distributions. The 
evolution of xmi is very similar in all models by construc- 
tion, while He n reionization is completed progressively later 
in models £1.6-a3 and £1.2-cv3, which have softer ionizing 
spectra. For more details on the results we refer the reader 
to C2012. 

It should be noticed that the outputs of the simulations 
are sufficiently close to each other in redshift space that 
there is always some overlap in the redshift range spanned 
by subsequent boxes. In particular, we will be using boxes 
Bi (with i = 1...54) which span the redshift range z ~ 6 — 15. 

While the above simulations are useful to investigate the 
effect of different reionization histories on the observability 
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Figure 2. Upper panel: evolution of the volume-averaged 21 cm 
optical depth for model £1.2-al.8 (black solid line), £1.2-al- 
3 (red dotted), £1.6-ec3 (blue dashed), £1.2-a3 (green dotted- 
dashed) and C4.39 (orang e long-dashed). Triangles refer to the 
values in the simulation bv lXu et al.l {2009). Bottom panel: evolu- 
tion of the volume averaged 21 cm optical depth in model £X.2-a3 
for different assumptions about the spin temperature, T s , being: 
the same as the gas temperature, T, determined from the radia- 
tive transfer simulation (green dotted-dashed line); the maximum 
between T and the temperature calculated including also the con- 
tribution from Lya photons from stars and x-ray photons from 
quasars (red dotted line). See text for more details. 



of the 21 cm absorption features, we have run an additional, 
higher resolution simulation in a 4.39/i _1 Mpc box with the 
aim of exploring instead the impact of resolution on our re- 
sults. This simulation has been run with 2 x 256 3 gas and 
dark matter particles, corresponding to a mass of the gas 
particle of 6.48 x 10 4 /i -1 Mq, and thus it is better suited 
to capture the absorption features due to high density gas 
at smaller scales. The simulation is similar to the one anal- 
ysed and discussed in Jeeson Daniel et al. (in prep.), but 
here we have chosen an amplitude of the total comoving hy- 
drogen ionizing emissivity £ = 1.8 (see eq. 3 in C2012) and 
a = 3 (in the remaining we will refer to this simulation as 
£4.39). These choices were motivated by the requirement 
that the evolution of the neutral fraction was similar to that 
of our reference run £1.2-a3. As a reference, simulation £1.2- 
q3 (£4.39) has xhi=0.971 (0.974), 0.519 (0.494) and 0.148 
(0.151) at 2=14, 9 and 7, respectively. Despite the fact that 
such a small box is not able to capture the reionization pro- 
cess as a whole, it still matches the available observables as 
the larger boxes do, and it is more appropriate to investi- 



gate resolution effects. Also in this case the outputs of the 
simulations are sufficiently close to each other in redshift 
space that there is always some overlap between subsequent 
boxes. Differently from earlier though, we have now a larger 
number of boxes Bi (with i = 1...450) spanning the same 
redshift range z ~ 6 — 15. 



3.1 Optical depth to the 21 cm line 

In the upper panel of Figure [2] we show the evolution of 
the volume-averaged (over the entire box) 21 cm optical 
depth for all models. The value of T2i cm is the same in all 
models at the highest redshifts, when the evolution of both 
the HI fraction and temperature are similar (see C2012 for 
more details). Differences are more evident at lower z. In 
particular, in model £1.2-a3 T2icm is up to one order of 
magnitude higher, because reionization is delayed compared 
to the other cases, for which differences are smaller. More 
specifically, £1.6-a3 has a similar xm, but lower temper- 
ature compared to £1.2-q1.8 and £1.2-q1-3, resulting in a 
slightly higher value of the optical depth. On the other hand, 
these latter two models have very similar values for both xm 
and T, and the differences in the volume averaged optical 
depth arise from more subtle details of the distribution of 
the physical quantities. As mentioned above, the reioniza- 
tion history in run £4.39 is very similar to that of £1.2-cv3, 
with a slightly (few percent) lower value of xm and a slightly 
(fewpercent) larger value of the volume averaged tempera- 



turqj. The above results in a lower T2i cm . 

In the lower panel of Figure [2] we show the evolution of 
the volume averaged 21 cm optical depth in model £1.2-a3 
for different assumptions about the spin temperature, T s . 
The green dotted-dashed line is our reference case, in which 
T s is assumed to be in equilibrium with the gas temper- 
ature, T, calculated from the radiative transfer simulation 
(which includes only UV photons). The same result is ob- 
tained when T s is calculated as the maximum between the 
gas temperature as determined from the radiative transfer 
calculation and that including also the contribution from 
Lya photons from stars. If in addition to the Lya photons 
we consider also the contribution from x-ray photons from 
quasars (red dotted line) , the gas temperature increases and 
as a consequence the optical depth decreases. The temper- 
ature with Lya and x-ray photons is calculated for a gas 
at the mean density using the semi-analytic approach as 
described in I Ciardi et al.l ( 2010T ). Here, the authors calcu- 
late the photon production from a population of stars and 
micro-q uasars in a set of cosm ological hydrodynamic simu- 
lations (jPelupessv et al. hoot ) which self-consistently follow 
the dark matter dynamics, radiative processes as well as star 
formation, black hole growth and associated feedback pro- 
cesses. We refer the reader to the original papers for more 
details. More specifically, the gas temperature as determined 
by Lya photons from stars is given by the short-dash-dotted 
blue line in Figure 2 of that paper, while the one with also 
heating from x-rays from quasars corresponds to the solid 



4 Despite having the same mean density of the larger box, the 
4.39/i — 1 box has gas at higher density. Because of this, more ion- 
izing photons are needed to maintain a similar ionisation fraction, 
with the consequence of an increased photo-heating. 
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Figure 3. Distribution of the 21 cm optical depth versus gas 
number density for model £1.2-a3 (left panel) and £4.39 (right). 
The rows refer to redshift z = 14, 9 and 7 (from top to bottom). 
The vertical lines refer to the volume-averaged gas number density 
in the box, which corresponds to S = 1. 



black line in the same plot. While in cells which have been 
(also partially) ionized the temperature is ~ 10 4 K, cells 
which have not been reached by ionizing photons have a 
temperature which can be as cold as that of the CMB. In 
this case, a contribution to the heating from Lya or x-ray 
photons can raise the temperature above its original value, 
reducing the optical depth to 21 cm. It should be noted 
that the above results depend on the model adopted for 
the Lya and x-ray photon production, and that this is not 
self-consistent with the RT calculations. While we will ad- 
dress this point more thoroughly in future developments 
by running radiative transfer simulations which follow self- 
consistently also the propagation of x-ray and Lya photons, 
the qualitative results will still hold. 

As a reference, the T2i cm calculated in model £1.2-a3 is 
~ 16% (50%) lower than the one found bv lXu et all ||2009h 
at z — 13.5 (9.3), and it becomes ~ 100% larger at z = 6. 
These values are lower than those found bv lMack fc Wvithd 
(2011) in their model without x-ray heating, most probably 
because they obtain an average gas temperature which is 
never above ~ 1000 K (also when x-ray heating is included) 
in the relevant redshift range, i.e. about one order of mag- 
nitude lower than the one found in our models. 

It should be noted that in all cases T2i cm is just an av- 



erage value, but the scatter is very large, as can be seen 
in Figure [3] where the distribution of 21 cm optical depth 
versus gas number density is shown for two reference simu- 
lations at redshifts z = 14 (top panels), 9 (middle panels) 
and 7 (lower panels). In the 35ft - Mpc box, most of the 
gas is initially concentrated either in high density and high 
optical depth cells (corresponding to cold, neutral, high den- 
sity regions) or in cells with T2i C m ~ lO^ 9 (corresponding to 
gas which has been ionized and heated to high tempera- 
tures). As reionization proceeds and more cells get ionized, 
the corresponding optical depth decreases. We also observe 
a larger spread of r2i cm with the gas number density, due to 
the fact that photons have now reached low density regions 
far away from the location of the sources of radiation. While 
at z = 9 there is still a large range of values of the optical 
depth, at z — 7 most cells have T2i cm ~ 10~ 10 , because 
reionization is almost complete. It is interesting to compare 
the distribution of the optical depth to the one of the gas 
temperature shown in Figure 7 of C2012. In fact, the two 
quantities exhibit a complementary behaviour, as expected. 
The distribution and evolution of T2i cm for £4.39 is similar, 
although the optical depth is always higher, and the range 
of densities spanned is larger compared to the other models. 

A more quantitative representation is given in Figure U 
where we show the percentage of cells as a function of the 
transmissivity, exp(-T2i C m), i.e. the quantity most relevant to 
observations (see following Sec). In all models, at the high- 
est redshift most of the cells have a high optical depth (i.e. 
exp(-T2icm) < 1), although, as already mentioned, there is a 
large spread and the cells which have already been ionized 
are concentrated around T2icm ~ 10 -9 (i.e. exp(-T2i C m) ~ 1). 
As ionization proceeds, the number of cells with such opti- 
cal depth increases, while those with higher T2i C m decreases. 
This results in a shift towards larger values of the transmis- 
sivity. At all redshifts there are also intermediate values of 
T"2icm, reflecting the distribution of densities and tempera- 
tures (see C2012). The limitations due to resolution effects 
are visible at all redshifts as a sudden drop of the number 
of cells below a minimum transmissivity. The behaviour is 
very similar in all models, with a systematically lower trans- 
missivity for £4.39. Due to the higher resolution of £4.39, 
larger values for the gas overdensity are sampled in this sim- 
ulation. It should be noted though that also in this smaller 
simulation the resolution is not enough to capture the even 
higher overdensities which are associated with collapsed ob- 
jects. This topic is the focus of a future investigation. 

In Figure [5] the gas overdensity of cells with optical 
depth larger than T2i C m is plotted at redshift z = 14, 9 and 
7 for models £4.39 and £T.2-a3. The shaded areas refer to 
the full range of overdensities, while the solid lines indicate 
the mean values. It is evident that, while the mean value 
is similar in both cases (because the hydrodynamical simu- 
lations have the same mean density), £4.39 reaches higher 
overdensities which produce the larger values of the optical 
depth. This is important for the observability of the line, as 
we will see in the following. 



4 MOCK SPECTRA 

From each simulation box Bi we extract 2500 random LOS. 
A total LOS is obtained as follows. We choose randomly a 
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Figure 4. Percentage of cells as a function of the transmissivity, 
exp(-T2icm), at 2=14 (upper panel), 9 (middle) and 7 (lower). 
The black solid, red dotted, blue dashed, green dotted-dashed 
and orange long-dashed lines indicate results for models £1.2- 
al.8, £1.2-al-3, £1.6-a3, £1.2-a3 and £4.39, respectively. 
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Figure 5. Gas overdensity of cells with optical depth larger than 
f"2icm at redshift z = 14 (upper panel), 9 (middle) and 7 (lower) 
for models £4.39 (orange shaded area in the back) and £1.2-a3 
(green shaded area in the front). The shaded areas refer to the 
full range of overdensities, while the solid lines indicate the mean 
values. 



LOS (Zosi) among those built from box Bi and we stack it 
with one (los?) chosen randomly among those built from box 
i?2. Whenever los\ and I0S2 overlap in redshift, we average 
the relevant quantities. For example, nm{z) of the total LOS 
is given by [nm,i(z) + riHi,i+i(z)]/2, where nm,i(z) is the 
Hi number density at redshift z in losi. This procedure is 
repeated for all the boxes. 

An example of such a LOS from model £ 1.2-a3 is shown 
in Figure [6] Here the panels (from top to bottom) indicate 
the evolution of H number density, Hi fraction, gas tem- 
perature and transmissivity, e _T21cm . The dotted lines indi- 
cate the volume-averaged values of the different quantities. 
It should be noted that here no smoothing has been ap- 
plied and the resolution shown corresponds to the one of 
the simulation, i.e. ~15-20 kHz depending on redshift. We 
will discuss the higher resolution simulations in Section [5J 
while here (Fig. [7]) we just show the transmissivity along a 
random LOS as a comparison with the larger box. Below, 
we will discuss the dependence of our results on a number 
of quantities, using the simulations run in the 35/i _1 Mpc 
box. 



4.1 Dependence on spin temperature 

In Figure [8] the dependence on different assumptions about 
the gas spin temperature is shown. The LOS is the same as 
that plotted in Figure[6] but now we show not only the refer- 
ence case (upper panel), but also a case in which T 3 is taken 
to be the maximum of the gas temperature as determined 
from the radiative transfer calculation (which includes only 
UV photons) and that including also the contribution from 
Lya photons from stars (middle panel), and from Lya plus 
x-ray photons from quasars (lower panel). More details on 
the calculation of the temperature are given in Section f3. II 
Because Lya heating is not extremely efficient, the resulting 
transmissivity is very similar to the case without such heat- 
ing. On the contrary, the contribution from x-rays is more 
substantial and this can clearly be seen in the Figure as a 
reduction of the transmissivity at redshifts when the extra 
heating is significant (z 12). This suggests that in the pres- 
ence of x-ray heating, the detection of the 2 1 cm forest would 
most pro bably not be fe asible (but see e.g. iMack fc Wvithel 
l201ll and lXu et al.ll201ll for a discussion of the dependence of 
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Figure 6. Physical state of the gas along a random LOS (solid 
lines) in model £1.2-a3. From top to bottom: physical H number 
density, H I fraction, gas temperature, gas transmissivity. The red 
dotted lines indicate the volume-averaged values of the different 
quantities. 
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Figure 8. Evolution of the gas transmissivity along a random 
LOS for different assumptions about the gas temperature, T. 
From top to bottom the gas temperature is: the one determined 
from the radiative transfer of UV photons only; the maximum 
of the one determined from the radiative transfer of UV photons 
and the one calculated including also the contribution from Lya 
photons from stars; the maximum of the one determined from the 
radiative transfer of UV photons and the one calculated includ- 
ing also the contribution from Lya photons from stars and x-ray 
photons from quasars. 
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Figure 7. Gas transmissivity along a random LOS in model 
£4.39. 



absorption features on the strength of the x-ray heating). It 
should be noted that the value of the x-ray background used 
here is appropriate for a population of high-z quasars. As a 
reference, at z = 10 it is ~ 4x 10 ~ 26 ergs cm -2 s _1 Hz -1 sr _1 
(see Fig. 1 of lCiardi et al]|201Ch . If the background were in- 
stead dominate d by a different class of sources, as x- ray bi- 
naries (see e.g. iFurlanettdkoOfj : iDiikstra et afll2012l ). then 
the associated heating would be different. We defer to future 
work a more detailed investigation of the role of the x-ray 
background and the contribution to it from different types 
of sources. 

As already poi nted out by other authors (e.g. 
iMack fc Wvithll201ll ). the above arguments seem to sug- 
gest that it would be possible to discriminate between dif- 
ferent IGM reheating histories, in particular if a high energy 
component were present in the ionising spectrum. 
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Figure 9. From the top to the bottom panels we show: the evolu- 
tion of the gas transmissivity along a random LOS for reionization 
model £T.2-al.8; the ratios between the transmissivity of model 
£1.2-al.8 and models £1.2-al-3, £ 1.6-a3, and £1.2-a3, respec- 
tively. The ratios have been calculated using the same LOS as 
the upper panel. 
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4.2 Dependence on source spectra 

As discussed in Section [3j the simulations have been run 
for a variety of source spectra. Here we show the differences 
induced in terms of transmission. In the top panel of Figure[9] 
a reference LOS for reionization model £1.2-al.8 is shown. 
In the other panels the ratios between the transmissivity of 
model £1.2-<ad.8 and models £1.2-al-3, £1.6-a3, and £1.2- 
q3 are plotted. In the third and fourth panels from the top 
the ratio is always larger than 1, meaning that the optical 
depth in £1.2-al.8 is lower than in £1.6-a3 and £1.2-a3 
(i.e. the simulations with softer spectra), as expected. In 
fact, although xm is very similar in £1.2-al.8 and £1.6-a3, 
the temperature in the latter is typically lower, resulting 
in a higher optical depth. On the other hand, as in model 
£1.2-a3 reionization is less advanced, the optical depth is 
always larger than for the other cases. Models £1.2-al.8 and 
£1.2-al-3 have very similar average xm and T, but their 
distributions are slightly different. This reflects also on the 
transmissivity. 

These very small differences in transmissivity predicted 
by different reionization models are to be expected. As dis- 
cussed in the previous Section, the models do have different 
values of T2i cm , but these are more pronounced when T2i cm 
is in any case too low to produce appreciable absorption. In 
fact, all models have very similar hydrogen reionization his- 
tories, with the largest discrepancies being associated with 
the temperature reached either in the ionized regions or in 
their outskirts, where the hardness of the spectrum makes 
most of the difference. From the point of view of the absorp- 
tion, though, the exact value of the temperature in a fully or 
partially ionized cell is not relevant. On the other hand, the 
regions which produce most of the absorption, i.e. the cold, 
high density gas not reached by ionizing photons, are basi- 
cally the same in all models. Thus, unless the reionization 
histories are very different from each other (e.g. for example 
a predominance of UV vs. x-ray sources), we do not expect 
to be able to distinguish them by means of observations of 
the 21 cm forest. 



5 INSTRUMENTAL EFFECTS 

In this Section we describe the process used to simulate the 
spectra as they would be observed by LDFAR. The input of 
the simulation is a file containing the optical depth as a 
function of redshift (or, equivalently, frequency) and the pa- 
rameters associated with the background radio source, i.e. 
its redshift, z s , flux density, Si n {z s ), an d spectral index, a . 
Consistently with previous studies (e.g. ICarilli et af1l2002T ) 
we assume that the source has a power-law spectrum. In 
principle, the effect of the Point Spread Function (PSF) side 
lobes running through the source of interest can be taken 
into account by including more sources at different posi- 
tions, with or without absorption features. Here, this has 
not been considered. For a narrow bandwidth, the equation 
giving the observed visibilities is: 



K(u) 



E 



r / \ — 27rzu-s , 

l v (s)e +7i, 



(2) 



(l,m,n) is a vector representing the direction cosines for a 
given source direction and n represents additive noise. The 
noise is given by the radiometer equation: 



1 



SEFD 



(3) 



where n s is the system efficiency, Av is the bandwidth, tint 
is the integration time and N is the number of stations. We 
assume that n s = 0.5 and N = A&. The system equivalent 
density is given by: 



SEFD 



2 K £} T S y S 

N dip ri a A c ff '' 



(4) 



min^f, 1.5626) 



where u = (u,v,w) are the coordinates of a given baseline 
at a certain time t, I v is the observed source intensity, s = 



where kjj is Boltzmann's constant, A e g 
is the effective area of each dipole in the dense and sparse 
array regimes respectively, iVdip is the number of dipoles per 
station (24 tiles times 16 dipoles per tile for a LOFAR core 
station) and r\ a is the dipole efficiency which we assume to 
be 1. The system noise T sys has two contributions: (i) from 
the electronics and (ii) from the sky. We assume that the 
sky has a spectral index of -2.55, obtaining T sys = [140 + 
60(^/150 MHz)" 2 ' 55 ] K. 

The complete Fourier plane sampling can be done by 
evaluating the above equation for every set of baseline co- 
ordinates. The predicted visibilities are then gridded and 
transformed via inverse Fourier transforms in order to ob- 
tain the dirty images. 

For the purpose of the 21 cm forest, fine spectral reso- 
lution is needed, which means that the spectra need to be 
predicted for a large number of channels (~ 10 000). Thank- 
fully, the computation can be done independently for each 
channel, and the maps can thus be evaluated in parallel. 
After assembling the full image cube, the LOS spectrum is 
extracted. 

In the upper panels of Figure [10] we show the spectrum 
of a source positioned at z s — 10 (i.e. v ~ 129 MHz). The 
intrinsic radio source sp ectrum, Sj n , is assumed to be simi- 
lar to Cygnus A (see also lCarilli et al.ll2002l : lMack fc Wvithd 
2011), with a power-law index a = 1.05 and a flux density 
Sin(zs) = 50 mjy. It should be noted that if the intrinsic 
radio spectrum were not a pure power-law, the identifica- 
tion of absorption features would be rendered even more 
complex by the requirement that the continuum level must 
be fit as well. The simulated absorption spectrum, Sabs, is 
calculated from the IGM optical depth found in our refer- 
ence simulation £1.2-q3. The observed spectrum, Sobs, is 
calculated assuming an observation time ti n t = 1000 h and 
a bandwidth Av = 20 kHz. The spectra are plotted at the 
resolution of the simulation, unless stated otherwise. While 
some absorption features are evident, it is not obvious that 
we will be able to disentangle them from the instrumental 
noise. Even though our simulation box reaching gas over- 
densities 5 as high as a few tens, typically a LOS through 
the IGM samples mild overdensities 8 ~ 1. As an example, 
while ^40% of cells along the LOS in the Figure (in the 
full redshift range z = 6.5 — 15) have T2i cm > 0.005, none 
has T2icm > 0.025, setting a limit on the g as transmissivity 
of e~y 21cm > 0.975. Thus, similarly to e.g. iMack fc Wvithd 
l|201ll ), we do not find narrow strong absorption lines (see 
below for further discussion). If we do a Gaussian smooth- 
ing of the spectrum over a scale s = 10 kHz (upper middle 
panel) or reduce the noise by a factor of 0.1 (similar to what 
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is expected from the SKA; upper right panel), then the ob- 
servability improves. This is more evident in the lower pan- 
els of the Figure, which show the quantity er a bs/o"obs, where 
cri = Si — Si u and i=abs, obs. 

Stronger absorption features appear if we use model 
£4.39, in which larger overdensities are present (see Sec- 
tion [3TTJ) . This is visible in Figure [TT] where Sm, Sabs and 
Sobs are shown for the same source of Figure 1101 The fre- 
quency range plotted is smaller to allow a better visualiza- 
tion of the absorption features. In this case the prospect for 
a detection improves, even halving the observation time (or, 
equivalently, the bandwidth). The third (fourth) panel from 
the left shows the case for Av — 2.5 kHz (Av = 10 kHz) 
and a smoothing of the spectrum over the same scale. While 
retaining a high resolution might result in a null detection, 
observation features are clearly visible if a larger bandwidth 
and smoothing scale are used. We have verified that features 
are detected for Av £ 5 kHz. Similarly to model £1.2-a3, 
along a LOS in the same redshift range, here we have ~40% 
of cells with r2i cm > 0.005, but we also find a handful of 
cells with T2icm > 0.1, corresponding to a gas transmissivity 
e -r 2 i cm < 0.905. 

If we were lucky enough to intercept one of such high 
density pockets of gas (with T2i cm > 0.1; these cells are 
found in ~ 0.1% of the LOS), very strong absorption fea- 
tures could then be detected, as shown in Figure [12] where 
we have located at z — 7.6 a radio source with the same 
characteristics as above. To calculate the observed spectrum 
we have now used a Av = 5 kHz and an integration time 
tint = 1000 h. An even better prospect of detection would 
arise if e.g. a high redshift Damped Lya System (DLA) were 
intercepted, but the simulations employed in this work do 
not have the resolution to investigate this possibility, which 
will be addressed elsewhere. 

Finally, moving towards higher redshift (z > 12), when 
most of the gas in the IGM is still neutral and relatively 
cold, would offer the chance of detecting a stronger average 
absorption. Although this range of redshift falls outside the 
one observed by LOFAR because of strong RFI contamination 
(Offringa et al. 2012, in prep), it is still interesting to look 
at the example shown in Figure 1131 where the usual refer- 
ence radio source has been located at z = 14. The Figure 
clearly shows that a LDFAR-type telescope could easily detect 
the global absorption (rather than the single absorption fea- 
tures observed at lower redshift) due to the highly neutral 
IGM if a relatively large bandwidth is used. In the Figure 
the results for Av = 20 kHz are shown, but a suppression 
of the source flux could also be detected for Av = 10 kHz. 
It should be noted, however, that distinguishing a suppres- 
sion of the source flux due to intervening neutral hydrogen 
from a source with an intrinsically lower flux would not be 
straightforward. 



6 CONCLUSIONS 

In this paper we have discussed the feasibility of the detec- 
tion of the 21 cm forest in the diffuse IGM with the radio 
telescope LOFAR. The optical depth to the 21 cm line has 
be en derived us i ng th e simulations of reionization presented 
m ICiardi et al.l (|2012h . which reproduce a number of obser- 
vational constraints such as the Thomson scattering opti- 
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Figure 12. Upper panel: Spectrum of a source positioned at 
z = 7.6 (i.e. v ~ 165 MHz), with an index of the power-law 
a = 1.05 and a flux density Si n (z s ) = 50 mjy. The lines are 
the same as those in Figure [lOl Here we have assumed the noise 
a n given in eq. [3] a bandwidth Av = 5 kHz, smoothing over 
a scale s = 5 kHz, and an integration time ti n t = 1000 h. The 
IGM absorption is calculated from the reference simulation £4.39. 
Lower panel: er a b s / cr obs corresponding to the upper panel. Note 
that this LOS has been chosen to intercept a strong absorption 
feature, with T2i cm = 0.12. 



cal depth dKomatsu et al an d the H i photo- ionisatio n 

rate at z ~ 6 (|Calverlev et alj20ld ; IWvithe fc Boltonlboill ). 
The simulations provide the evolution of the spatial distri- 
bution of relevant physical properties such as gas neutral 
fraction, temperature and velocity field. The main results of 
our investigation can be summarized as follows. 

• The spectra from reionization models with similar total 
comoving hydrogen ionizing emissivities, but different fre- 
quency distributions look remarkably similar. Thus, unless 
the reionization histories are very different from each other 
(e.g. a predominance of UV vs. x-ray heating), we do not 
expect to distinguish them by means of observations of the 
21 cm forest. 

• The photo-heating associated with the presence of a 
strong x-ray background would make the detection of the 
21 cm line absorption impossible. The lack of absorption 
could then be used as a probe of the presence/intensity of the 
x-ray background and the thermal history of the universe. 

• Along a random line of sight LOFAR could detect a global 
suppression of the spectrum from z Si 12, when the IGM is 
still mostly neutral and cold, in contrast with the more well- 
defined, albeit broad, absorption features visible at lower 
redshift. Sharp, strong absorption features associated with 
rare, high density pockets of gas could be detected at z ~ 7 
along preferential lines of sight. 
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Figure 10. Upper panels: Spectrum of a source positioned at z s = 10 (i.e. u ~ 129 MHz), with a power-law index of a = 1.05 and a 
flux density Si n (z s ) = 50 mjy. The red dotted lines refer to the intrinsic spectrum of the radio source, Si n ; the blue dashed lines to the 
simulated spectrum for 21 cm absorption, S a b s ; and the black solid lines to the spectrum for 21 cm absorption as it would be seen by 
L0FAR. The latter has been calculated assuming an observation time tint = 1000 h and a bandwidth Au = 20 kHz. The IGM absorption 
is calculated from the reference simulation £1.2-03. Left panel: S a b s and Sobs are obtained without smoothing, s = 0, and with the noise 
tin given in eq. \3\ Middle panel: S a bs and Sobs are obtained after smoothing over a scale s = 10 kHz, and with the noise a n given in 
eq. \3\ Right panel: S a b s and S b s are obtained without smoothing, s = 0, and with l/10th of the noise a n given in eq. [3] Lower panels: 
"abs/^obs corresponding to the upper panels. Here Oi = Si — Si n , with i = abs, obs. See text for further details. 
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Figure 11. Upper panels: Spectrum of a source positioned at z = 10 (i.e. v ~ 129 MHz), with an index of the power-law a = 1.05 and 
a flux density 5i n (2 s ) = 50 mjy. The lines are the same as those in Figure [lOl S D b s has been calculated assuming the noise a n given in 
eq. [3] The IGM absorption is calculated from the reference simulation £4.39. From left to right the panels refer to a case in which S a b s 
and Sobs are obtained: with a bandwidth Av = 20 kHz, without smoothing, s = 0, and with an integration time ti n t = 1000 h; with 
a bandwidth = 20 kHz, without smoothing, s = 0, and with an integration time tj nt = 500 h; with a bandwidth Au = 2.5 kHz, 
smoothing over a scale s = 2.5 kHz, and with an integration time tint = 1000 h; with a bandwidth Au = 10 kHz, smoothing over a scale 
s = 10 kHz, and with an integration time ti n t = 1000 h. Lower panels: o" a b s /a- b s corresponding to the upper panels. 
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Figure 13. Upper panel: Spectrum of a source positioned at 
2 = 14 (i.e. v ~ 95 MHz), with an index of the power-law 
a = 1.05 and a flux density Si n (z s ) = 50 mjy. The lines are 
the same as those in Figure llOl Here we have assumed the noise 
a n given in eq. [3] a bandwidth Au = 20 kHz, smoothing over 
a scale s = 20 kHz, and an integration time ti n t = 1000 h. The 
IGM absorption is calculated from the reference simulation £4.39. 
Lower panel: <T a b s /<Tobs corresponding to the upper panel. 



The most challenging aspect of the detection of a 
21 cm forest remains the existence of high-z radio loud 
sources. Although a QSO has been detected at z = 7.085 
|Mortlock et al.ll201ll ). the existence of even higher redshift 
quasars is uncertain. The predicted number of radio sources 
which can be used for 21 cm forest studies in the whole sky 
per unit redshift at z = 10 varies in the range 10 - 10 4 de- 
pending on the model adopted for the luminosity function 
of such sources and the instrum ental characteristics (e.g. 
ICarilli et al.ll2002l ; IxiT et al. 2009), making such a detection 
an extremely challenging task. Th e possibility of using GRB 
afterglows has been suggested by lloka fe Meszaros! J.2005), 
concluding that it will be difficult to observe an absorption 
line, even with the SKA, except for very energetic sources, 
such as GRBs from the first stars. In fac t, a similar calcula - 
tion has been repeated more recently bv lToma et ail [|20ld ) 
for massive metal-free stars, finding that the flux at the same 
frequencies should typically be at least an order of magni- 
tude higher than for a standard GRB. 

An absorption feature stronger than the one produced 
by the diffuse IGM, would be the on e due to intervening star- 
less m inihalos or dwarf galaxies (i.e. IXu et al.ll201ll ; lMeiksinl 
l201ll ), resulting in an easier detection. On the other hand 
the optical depth would strongly depend on the feedback 
effects acting on such objects. Because of the large uncer- 
tainties about the nat ure and intensity of hi gh-z feedback 
effects (for a review see lCiardi fc Ferraral 2005 and its ArXiv 
updated version), it is not straightforward to estimate the 
relative importance of these two absorption components un- 



less a self-consistent calculation is performed. We defer this 
investigation to a future paper. 
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